Flux penetration in slab shaped Type-I superconductors 
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We study the problem of flux penetration into type-I superconductors with a high demagnetization 
factor (slab geometry). Assuming that the interface between the normal and superconducting regions 
is sharp, that flux diffuses rapidly in the normal regions, and that thermal effects are negligible, we 
analyze the process by which flux invades the sample as the applied field is increased slowly from 
zero. We find that flux does not penetrate gradually. Rather there is an instability in the process 
and the flux penetrates from the boundary in a series of bursts, accompanied by the formation of 
isolated droplets of the normal phase, leading to a multiply connected flux domain structure similar 
to that seen in experiments. 
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When a type I superconductor is placed in a magnetic 
field less than the bulk upper critical field it exhibits a 
phase with interpenetrating domains of the normal and 
superconducting state called the intermediate state. De- 
spite over fifty years of work, a complete description of 
the physics of this phenomenon has been elusive. The 
earliest theoretical work on this problem goes back to 
Landau jlj who studied the problem of the equilibrium 
configuration of an infinite superconducting slab of thick- 
ness d placed in an applied magnetic field B app . When 
d is much larger than the penetration depth A and the 
coherence length £, he showed that the Meissner state is 
unstable to a configuration of alternating superconduct- 
ing and normal laminae. Other configurations of flux do- 
mains, various sample geometries, etc. have been studied 
both experimentally and theoretically |jj(and refer- 
ences therein). 

It has been clear since the work of Landau that these 
structures arise from a competition between the mag- 
netic field energy, the condensation energy of the super- 
conducting regions, and the surface energy of the inter- 
face between the normal and superconducting regions. 
However, the striking fact about most of the experiments 
is that the regular structures envisioned by Landau are 
rarely seen. Instead one sees complicated patterns which 
are strongly dependent on the temperature, disorder and 
field history of the sample. 

In the last few years considerable progress has been 
made in understanding a related problem: flux pene- 
tration into (or expulsion from) long cylindrical samples 
oriented parallel to the magnetic field Q . Although in- 
teresting flux patterns are seen in these samples as tran- 
sients in going from a normal to a superconducting state, 
a steady-state intermediate state, which arises due to de- 
magnetization effects, is not seen. Very recently, Gold- 
stein, Jackson and Dorsey (hereafter referred to as GJD) 
Q attempted to analyze the influence of demagnetization 
on the formation of flux structures for slab-like type-I 



superconductors. They assumed that transient currents 
in the normal domains decay very fast (which amounts 
to the normal state conductivity being zero) , and argued 
that one can then think of the dynamics as a simple gra- 
dient descent of the free-energy. When the coherence 
length £ is small, the interface between normal and su- 
perconducting regions is approximately sharp. The free- 
energy for a sample placed in an applied magnetic field 
B ap p can then be written as 

T(B app , A) = F B +F C + T s . (1) 

Here the first term denotes the magnetic free en- 
ergy, the second term the condensation energy which is 
—H^/8tt per unit volume in the superconducting regions, 
and the third the interfacial energy which is H^AA/8ir. 
A is the surface energy parameter which is of the order of 
the coherence length £, H c is the bulk upper critical field 
and A is the area of the interfaces between the normal 
and superconducting regions. 

GJD p| assumed that (in addition to a bulk term) the 
magnetic part of this energy could be written as a long- 
range interaction between current loops localized on the 
interfaces between the normal and superconducting re- 
gions. This made the problem very similar to the problem 
of the dynamics of two-dimensional fcrrofluid droplets in 
a magnetic field || . In the ferrofluid case it is known || 
that regular shapes evolve continuously into labyrinthine 
patterns when the applied field is increased adiabatically. 
In the case of superconductors, GJD showed that a cir- 
cular flux droplet in a sea of superconducting material 
(with an area much larger than the equilibrium area) 
changes into a many armed structure with three-fold co- 
ordinated nodes. Such convoluted structures are indeed 
seen in some experiments 0] . They also calculated the 
equilibrium periodicity for a laminar structure and got re- 
sults that are numerically close to Landau's. While this 
would suggest that their free-energy captures the physics 
of this situation, GJD [^j were unable to derive it from 
the more basic Ginzburg-Landau description. 
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As shown below, a careful analysis of the magnetic 
free energy yields an interaction between superconduct- 
ing domains that cannot be expressed in terms of current 
carrying loops (at least for low fields) , and results in qual- 
itatively different behavior. (Late stages of flux invasion, 
when the supercondcting domains are tall and thin, are 
discussed later.) It is not difficult to see why supercon- 
ductors and ferrofluids are so different. When a thin layer 
of a ferrofluid of thickness d and a characteristic trans- 
verse linear extent of L is placed in an applied magnetic 
field B app , to leading order in d/L the B field inside (and 
outside) the ferrofluid is equal to B app . This uniform 
field induces a uniform magnetization within the sample, 
which (using the equation VxB = —jext + 4-7rV x AI) 
gives rise to a small fringing B field near the sample 
edges. (The fringing field in turn induces a higher order 
non-uniform magnetization.) To leading order in d/L 
the induced currents consist of a ribbon flowing around 
the boundary of the ferrofluid, causing a long-range in- 
teraction between different parts of the boundary. On the 
other hand, for a thin superconductor in the intermedi- 
ate phase, B is not even approximately equal to B app 
near the superconducting regions. Outside the sample, 
just above or below a superconducting domain, B is par- 
allel to the surface, while inside the domain B is zero. 
In addition to ribbons of current along the side walls of 
the superconducting regions, there are also large current 
sheets on the top and bottom surfaces, dominating the 
inter-domain interaction. Thus while the basic GJD idea 
of long-range interactions destabilizing regular patterns 
is correct, the actual description of the experimental pat- 
terns is more complicated. The numerics we report in this 
paper show flux invading in bursts, pinching off from the 
boundaries to form droplets. This is qualitatively differ- 
ent from the continuous evolution of GJD H . 




FIG. 1. B-field for a slab shaped (a) superconducting do- 
main (b) ferrofluid in a vertical applied field. The distortion 
of the field is small only for (b). 

We now consider Tb in detail. For an arbitrary sample 
placed in an applied magnetic field B app , the magnetic 
free-energy is T B = ^ J d 3 x(B 2 - 2B.H). Here B is 
the magnetic induction, H is the magnetic field and the 
integral is over all of three dimensional space. It is con- 



venient to cast this equation in a slightly different form: 
since VxH = V x B app , and B is a transverse field, H can 
be replaced by B app above. Adding the £?-independent 
term B 2 pp /8ir to the free energy density yields 

^b = -^J d 3 x(B - B app ) 2 (2) 

Thus in order to evaluate Tb, one has to find the mag- 
netic induction B for a given applied field B app . For 
the superconducting regions, where B — 0, the contri- 
bution to Tb is simple. The non-trivial part of the 
calculation consists of determining the magnetic field 
in the normal domains and in the space outside. Out- 
side the superconducting regions, one can define a mag- 
netic scalar potential <f> by B — B app = v</>> satisfying 
V 2 = 0. (We assume that the applied magnetic field 
is not changed rapidly, so that transient currents in the 
normal regions can be ignored.) Given any configuration 
of superconducting and normal regions we then have to 
solve Laplace's equation outside the superconducting re- 
gions. Since the normal component of B is zero at the 
boundaries, for a flat superconducting slab in a vertical 
field B app pointing upwards the boundary conditions are 
(i) d n cj) — ±B app on the top and bottom surfaces of the 
superconducting regions and (ii) d n <p = on the inter- 
faces between the normal and superconducting regions, 
where d n is the normal derivative. (We have assumed 
that the interfaces are vertical, ignoring fanning out of 
flux domains at the top and bottom surfaces.) 

At this stage, we use the quasi two-dimensional na- 
ture of the problem to reduce its computational com- 
plexity. For a thin sample, most of Tb, jp J d 3 x(\7(j)) 2 , 
is stored outside the sample (and in the normal re- 
gions). Integrating by parts transforms this expression 
into J s ds(f)d n (f>, where the integral runs over the top 
and bottom surfaces of the superconducting regions. (On 
the side walls, the normal derivative of cj> is zero.) In eval- 
uating the surface integral, for a thin sample the top and 
bottom surfaces of the superconducting regions can be 
treated as approximately coincident, lying in the z = 
plane. By symmetry, <j){z = 0) = 0, except in the super- 
conducting regions where <p(z = + ) = — <\>(z = CP) ^ 0. 
Using the boundary condition d z <fi(z = ) = —B app for 
the superconducting regions yields 

Fb = ^ jd^(z = + ,f) (3) 

where f is a two-dimensional vector. The integral is 
performed only over the superconducting regions. To ob- 
tain <f)(z — + ) we solve the inverse problem: if <p(z = + ) 
were known, d n (fi(z = + ) could be found by solving 
Laplace's equation with Dirichlet boundary conditions. 
Using d z (p(z = + ) = —B app yields the condition 
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Ba W = -0.M» = O + ,f)= / dr , K(r-7)4>(z = + ,r'). 

(4) 

The integral on the right hand side is over the supercon- 
ducting regions. K is given by 

1 2z 2 — \r — r'l 2 
K(r-n=^a- ¥ _^ 2 + z2]5/2 (5) 

whose (two dimensional) Fourier transform is given by 
|fc|/27r. Since K is convoluted with 0, which is a smooth 
function, it can be replaced by 

K(f-r') = -l/[27r|f-f | 3 ] \r-r'\>e (6) 

with a compensating (5-function of strength l/e at the 
origin, in the limit e — > 0. (For the numerics presented 
in this paper, the cutoff e is effectively the lattice size.) 

For an arbitrary pattern of superconducting and nor- 
mal regions, one can only evaluate Tb numerically using 
the prescription above, ft is not possible to reduce the 
expression to integrals over domain boundaries. (Note 
that it is not possible to solve Eq.(^) trivially by Fourier 
transforming, since the right hand side is only integrated 
over the superconducting regions.) The special case of a 
single circular superconducting region can be solved ana- 
lytically and can be used as a check of the numerics. 

We now discuss the numerics. Since the experimen- 
tal patterns depend a great deal on the field history of 
the sample, we concentrate on the following question: if 
one slowly increases the field from zero, how does flux 
penetrate the sample? We choose a sample with a thick- 
ness d ~ 12A and linear dimension L = lOd (as stated 
before A is essentially the coherence length £). In ex- 
periments the flux domains branch near the surfaces for 
d > O(800A), and type-II behavior is seen for d < 0(A). 
Our choice of d avoids both these regimes. Because 
L >> d, we can use the two dimensional formulation 
of the magnetic boundary value problem to find Tb ■ A 
typical value for type-I superconductors, A ~ 1500A, 
corresponds to L ~ 15/im. Although this is much smaller 
than experimental sample sizes, the qualitative aspects of 
our results should apply to larger samples as well. 

We divide the sample into lattices of various sizes; for a 
61 x 61 lattice the lattice constant is approximately 2.4A. 
The plaquettes in the lattice are either superconducting 
or normal. Apart from the isolated normal droplets, the 
superconducting-normal boundary consists of an outer 
interface close to the sample edges. There are two basic 
moves in the numerics: changing any one superconduct- 
ing plaquette on the outer interface to a normal plaque- 
tte, or vice versa. (The dynamics of the droplets are dif- 
ferent, due to flux conservation, and are discussed later.) 
Both possibilities need to be allowed for in order allow 
for droplet formation. A plaquette is flipped if there is 
a force that favors the move; the force is the gradient of 



the free energy, which is given by the sum in Eq. (|l|) . The 
surface tension force (from ,F s )is calculated by smoothing 
the lattice interface and computing the local curvature. 

Trying single plaquette moves requires a special treat- 
ment of surface tension in the dynamics, since flipping a 
single plaquette corresponds to a small sharp protrusion 
of the interface, with roughly the same forward and lat- 
eral extent. For a small lattice constant such a move will 
always cost a large surface energy compared to the en- 
ergy gained from the magnetic and condensation terms. 
Thus in the lattice dynamics whenever the force on a 
segment of the interface favors flipping a plaquette, it is 
likely that at the next time step it will be favorable to 
revert the plaquette to its original state. This is because 
in a lattice approximation the interface is forced to make 
larger excursions than it would like to; it would be better 
for a segment of the interface to move forward by only 
a fraction of a plaquette, and then let neighboring seg- 
ments catch up. We use a simple prescription to cure 
this lattice effect: when a plaquette is flipped, it is not 
allowed to flip back at the very next time step, although 
it can flip back thereafter. This should give rise to errors 
in the pattern only of the order of a lattice constant. 

We start with B app = 0, when the whole sample is su- 
perconducting. We raise the field till one plaquette on 
the boundary becomes normal. Flipping this plaquette 
can make it favorable to flip other plaquettes, in which 
case we let the system evolve till it reaches a stationary 
state. At this point the field is raised again. This pro- 
cess simulates the adiabatic increase of the magnetic field 
that we wish to study. Isolated droplets have to be han- 
dled differently, since the flux in them is conserved. In 
our simulations this constraint is obeyed approximately: 
once a droplet is formed, we move it rigidly in the direc- 
tion of the force on it till the force is zero. At this stage 
we adjust the number of plaquettes in the droplet so that 
the flux in it is as close as possible to its original value. 

We now discuss our results. When the field is raised 
just above the lower critical field, one plaquette becomes 
normal. If the penetration of flux were gradual, one 
would expect to have to raise the field further for more 
flux to enter. Instead we found that at a field just 
slightly above the lower critical field, the flux penetrates 
a distance into the sample of the order of twenty times 
the coherence length before the first droplet pinches off. 
Increasing the field further produces similar behaviour: 
much of the evolution is in the form of bursts of magnetic 
flux penetrating from the boundary which then pinch off 
to form droplets. This reflects an instability in the pro- 
cess of flux penetration and is the main result of our 
work. In Fig. |^ we show the patterns seen for L = 61 at 
a field only moderately above the field of first flux pene- 
tration. We saw similar patterns at comparable fields for 
the other lattice sizes. The droplets form near the bound- 
ary of the sample and then move towards the center of 
the sample, leaving a region near the boundary flux-free, 
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similar to what is seen in experiments. We also saw that 
the droplets typically shrink when one puts in flux con- 
servation, though this observation may not have much 
experimental significance in view of the fact that we deal 
with the collective motion of the droplet approximately. 




conducting regions become thinner at higher fields, it is 
not clear whether quasi two-dimensional descriptions are 
valid here. 

Although the existence of droplets is fairly ubiquitous 
in the experiments, the actual shape of the droplets varies 
rather widely, from compact droplets in experiments on 
mercury in the early stages of flux penetration Jfj], to 
long laminar structures in experiments on lead JlCQ . In 
order to explore this further, one would need to satisfy 
the constant flux constraint accurately, as well as con- 
struct a continuum description (presumably analytical) 
of the dynamics. A continuum analysis would also be 
necessary to treat surface tension exactly. 

To conclude, we have developed a description of the 
problem of flux penetration into slab shaped Type-I su- 
perconductors based on the sharp interface approxima- 
tion. Numerical simulations on a lattice show that as 
the applied magnetic field is increased, flux penetrates 
in bursts, forming droplets, leading to isolated normal 
regions. While the multiply connected nature of the pat- 
terns has been emphasized in the literature || , the insta- 
bility that we have noticed does not seem to have been 
reported so far. This instability should be apparent in 
real time imaging of the process of flux penetration. 
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FIG. 2. Droplet state for (a) numerical simulations on a 
61 x 61 lattice (b) experiments on mercury |J. Normal re- 
gions are black in (a) and white in (b). The flux front near 
the top in (b) is due to finite sample thickness |ll| 

The formation of droplets is promoted by the discon- 
tinuous nature of the flux invasion, which causes rela- 
tively large amounts of flux to enter the system at low 
magnetic fields. A large amount of flux penetration re- 
duces the magnetic forces which tend to drive further flux 
into the system. It is then possible for a normal region 
to find it favorable to revert to being superconducting, 
which is how droplets are formed. At later stages when 
the applied field is higher such reversion to a supercon- 
ducting state becomes less likely. Further flux coming 
in from the boundary would then probably coalesce with 
the already existing droplets, leading to the labyrinthine 
patterns seen in some experiments, reminiscent of the 
patterns obtained by GJD. || However, since the super- 
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